A DISCRETE INVERSE SCATTERING TRANSFORM FOR QS^ 
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Abstract. We derive a fully discrete Inverse Scattering Transform as a method for solving the initial- value 
problem for the Q3|5 lattice (difference-difference) equation for real-valued solutions. The initial condition is 
given on an infinite staircase within an N-dimensional lattice and must obey a given summability condition. 
The forward scattering problem is one-dimensional and the solution to QS^ is expressed through the solution 
of a singular integral equation. The solutions obtained depend on N discrete independent variables and N 
parameters. 



1. INTRODUCTION 

The Q3s equation 

—] 

iPQj 

where = {p^ — a?){p^ — b^) and = {q^ — a^){q^ — b^), is a nonhnear lattice (difference-difference) 
equation. The dependent variable u depends on two discrete independent variables n, m S Z in the following 
way 

u — u{n^m) M = u(n+l,m) M = M(n, m-|-l) u = u{n + \,m -\- 1). 

The lattice parameters p, q are associated with the n- and m-directions respectively, and a, 6, 5 are additional 
parameters. In 2003 Equation (jl.ip appeared in the ABS classification list j5| of such scalar affine-linear 
partial difference equations defined on four points, which possess the multidimensional consistency prop- 
erty [23 [7]. The particular parametrisation given here is due to [25] who also gave 7V-soliton solutions to 
Ijl.ip using a Cauchy matrix approach (see also [23] [H])- The extra parameter 5 may be set to 1 with- 
out loss of generality, and in the special case (5 = equation (|l.ip can be (gauge) transformed to the NQC 
equation [26j , which is a particular discretisation of a degenerate form of the Krichever-Novikov equation [20] . 



P{uu~\- uu) — Q{uu + uu) — {p^ 



{uu 



(1.1) 



Equation p.l|) defines a nonlinear relation between four points of the lattice shown in Figure (U and given 
its affine-lincarity, may be solved uniquely for each point in terms of the other three. 



X: 




m 



-> n 



Figure 1. Elementary lattice quadrilateral 

This property allows one to define a well-posed initial value problem [5] [25] [2S] by giving an initial profile 
along an infinite staircase within the (n, m)-lattice (see Figure l^al) . Due to the multidimensional consistency 
of (jl.ip however, this concept can be generalised. If we consider QSa living in a multidimensional lattice, 
with copies of the equation imposed on each elementary quadrilateral within this space, then we can still set 
up a well-defined initial value problem by giving an initial profile along some multidimensional staircase 
r. This is shown in Figure [2bl 
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(b) Multidimensional staircase F of 
(a) Staircase of initial values in the initial values 

(n, m)-lattice 

Figure 2 



In this paper we solve the inital value problem for (jl.ll) for an initial condition given along a multidimen- 
sional staircase F, within an N-dimensional lattice. 



1.1. History. 

The Inverse Scattering Transform (1ST) has been widely used as a mathematical tool for obtaining solutions 
of integrable nonlinear partial differential equations since its discovery in the 1960s by Gardner, Greene, 
Kruskal and Miura [16] [17] . Here the authors used it to solve the initial- value problem for the KdV equation, 
and the method has since seen many generalisations. Some other physically relevant partial differential 
equations which are also solvable by the 1ST include the nonlinear Schrodinger equation [33^ and the sine- 
Gordon equation [1]. 

The first application of the 1ST to differential-difference equations dates back to 1973 with Case and Kac 
[Hj [H], who considered a discretisation of the Schrodiner equation on the half-line n > 0. In 1974 Flaschka 
[14] showed how this could be applied to solutions of the Toda lattice and Ablowitz and Fadik [2] [3] [4] 
then derived a new dicrete scattering problem and showed how it was applicable to a number of physically 
important systems. 

More recently in 1999 Shabat [31] found a new discrete scattering problem by applying Darboux trans- 
formations to the Schrodinger equation. This problem has since been studied by various authors [5] [SO] [21] 
for the cases of continuous and discrete time evolution, and as an eigenvalue problem for analytic difference 
operators. This scattering problem also appears in the rigorous formulation of the discrete 1ST for the lattice 
potential KdV equation given in [10], and was also generalised to a multidimensional setting in [9]. While 
[9] does consider the 1ST for Q3s, there it was required that the "potential" term arising from the initial 
condition decay exponentially. Here however the potential term arising from the initial condition must only 
decay faster than n^^, which greatly widens the class of solutions obtained. 



1.2. Outline of Results. 

In this paper we solve the initial value problem for the Q3s equation (|l.ip by developing a fully discrete 
1ST for this equation which naturally incorporates its multidimensional consistency. We begin in Section [2| 
by giving the Lax pair for (|1.1|) and looking at the associated linear problem arising from this. In Section |3| 
we then set up the 1ST framework (reality assumptions, initial value space, boundary conditions, etc.) within 
the multidimensional lattice. The forward scattering problem is carried out in Section|4l where Jost solutions 
are constructed as functions existing on the staircase F. An important step is then given in Section [S] where 
we make assumptions on the boundary conditions for the solutions off F. This is essentially imposing that 
our boundary conditions are "time" independent. The inverse problem is then carried out in Sections [6| and 
\7\ We then show in Section [5| how to reconstruct the solution u of Q3s- A key part of this procedure is the 





( Pu 




y uu 




f Qu 


u 1 


\ uu 




A DISCRETE INVERSE SCATTERING TRANSFORM FOR Q3s 3 

solving of the singular integral equation 

for the function ^±fc(C), where the plane- wave factors are contained in p and all other quantities appearing 
in the equation are known in terms of the scattering data. Finally a one-soliton example is given in Section 
[9l For convenience some of the longer proofs are given in the Appendix. 

2. Linear Problem for Q3a- 
A Lax pair for equation (II. ip . which was first given in |24j . is 

(2.1a) 

^ ^ ^^^„ ^ , (2.1b) 

where this system is consistent (i.e. = ^) if and only if u solves (jl.ip . Here C is the spectral parameter. 
The dual function U is determined by solving the first-order equations 

UU = P{u^ + u^) - (2p2 - a2 _ 52^)^ ~ _^ _ (2.2a) 
UU = Q{u^ + u^) - (2g2 - - b^)uu+ —. (2.2b) 

These equations can be solved to determine U. The integration constants are determined by the boundary 
conditions that will be imposed on U. This is the setup for the linear problem of QS^- in the (n, m)-plane. 
In order to incorporate the multidimensional consistency of this equation however, we consider how the Lax 
equations ()2.ip change when we we make a shift in any one of the N possible directions within the lattice. Let 
the N discrete independent variables be denoted by ni, tt-n (all elements of Z), with associated parameters 
Pi,...,Pn- If we choose an arbitrary direction with independent variable and associated parameter p^, 
then we let u denote a shift of u in this rifc-direction, that is 

u = u(ni,n2, ...,nfc, ...,nN]Pi,P2, —,Pk, -Pn) 
u^ u{ni,n2, + 1, ...,nN\Pi,P2, -jPk, —Pn)- 

Due to the symmetry of the equation and its multidimensional consistency, the Lax equation is this n/j- 
direction is 

f p,u-ipi~b^)u e-h^ \ 

W ^ UU~±§^^ ipl-b^)u-P,u I 

Equation (|2.3p is then compatible with every other of the N-1 remaining Lax equations (we have ^N(N-l) 
Lax pairs), provided that u solves Q35 in each pair of lattice directions. The dependence of U on is found 
by solving 

UU = Pk{u^ + u^) - {2pl -c? ~ b^)uu+ —, (2.4) 
and again using the boundary conditions of U to determine the integration constants. 



3. Setting up the Discrete 1ST 
Here we set up the discrete 1ST in the multidimensional lattice. 
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3.1. Reality Assumptions. 

While equation (|l.ip is defined for complex- valued solutions u, for the purposes of the 1ST we assume that 

- The solution u and all parameters a, b,pi, ...,pn are real. 

- Since the equation only depends on the squares a^, , we choose to set a > 5 > and 
Pk>0 for aU k = 1,...,N. 

Complex- valued solutions were allowed in [9], however this led to a significantly stronger restriction on the 
initial condition. Our aim here is to obtain the widest possible class of solutions. 



3.2. Staircase. 



We now define the staircase F of initial conditions. Let us first consider some examples. 

Suppose first that we choose to give the initial condition on a line To spanned by the variable n (and 
parameter p), which corresponds to one of the N variables Uk (with parameter pk). This is shown in Figure 
[5al This is the usual setup for the continuous 1ST for equations such as the KdV equation. Along this line all 



Fi 



Fo 



(a) Line Fo of initial values 



(b) (1, 1) staircase Fi of initial val- 
ues 



Figure 3 



other N-1 variables are held constant. From the Lax equation p.3p . by eliminating the second component, 
the first component of <p satisfies 

The second component in then constructed from 

(C' - 02 = (j? - C')^ It ?1 - [PU- (/ - \?)U) (/.I. 

This gives as a function of n along the line Fq, which is the direct scattering problem. In order to obtain 
as a function of the remaining N-1 lattice variables, we need to give boundary conditions at rt — )• —oo 

01 n +0O, and then use the remaining N-1 Lax equations (|2.3|) to determine the "time" evolution of the 

scattering data in each of these directions. In the language of the continuous theory this discrete 1ST is a 

l+(N-l)-type scattering problem. 

Contrary to the continuous case, one of the benefits of the discrete setup is that we can easily change 

the one-dimensional manifold along which we specify the initial conditions. Perhaps the most natural such 

manifold is a (l,l)-staircase in the (n, m)-plane, shown in Figure l3bl For each 3-point segment of the staircase 

which iterates first in the n-direction and then in the m-direction we have 

(.^ - c^)^ I - (91^(1^}1^^ + (,^ _ . 0, 
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while for each 3-point segment which iterates first in the m-direction and then in the n-direction 



P u ~ {p^ - q^)u — Qu 



0. 



Solving these equations gives 01, and then 02 may be constructed from similarly considering the first com- 
ponent of the two Lax equations. Thus the result of the forward scattering problem is that we know 0) as a 
function along Fi, that is in terms of some independent staircase variable which depends on n and m. If we 
let this new variable be denoted by i, and let io correspond to the point (no, mo) on Fi, then assuming that 
we iterate first in the n-direction, the change of variables from i to n,m is given by 



1 , 

2(^ + 1 



«o) 



(3.1) 



where the brackets denote the floor function. This is once again a l-|-(N-l)-type scattering problem. 

We now look at how this may be generalised. Consider an arbitrary staircase F which has the following 
properties: 

- F is an infinite staircase which iterates in I of the N lattice directions, where 1 < Z < N 

- Every iteration along F corresponds to a positive iteration in one of the lattice variables n/j 

- F is defined through some repeated stepping algorithm. 

Figures l4al and l4bl show examples of such staircases, where we assume that both of these are repeated infinitely 
in both directions. Figure Hcl shows an example of a staircase which violates the above criteria. 




(a) (2, 2) staircase of initial 
values 



(b) Repeated three- 
dimensional staircase of 
initial values 



(c) Inadmissible staircase of 
initial values 



Figure 4 



3.3. Linear Problem. 

Given a staircase F of initial values, we now look at the linear problem along F. Since each iteration along 
F is in one of the lattice directions, the Lax equation for each iteration will be one of the I equations (|2.3p . 
From here it is convenient to introduce the staircase variable i, which will cycle through the various lattice 
variables encountered along F, and will be related to these I lattice variables by a relation similar to p.ip . 
To take care of the fact that the parameters pk will also change along the staircase we introduce the staircase 
parameter p = p(i) (and = (p^ — a^)(p^ — fo^)), which will cycle through the parameters pk encountered 
along the staircase. This new variable i and parameter p(i) allow the scattering problem for 0i (obtained by 
eliminating 02) along F to be conveniently expressed as 

(P^ - e)^^ - (^iiJi!:^!)^) + _ c^)^ 0, . 0, (3.2) 

where 

01 = C), 0i = 0i(« + i;C), p = p(i)7 p = p(« + i)- 

The second component of the eigenfunction is then constructed from 

(C' - &') h - (p' - C')^ 1101 - {Pu - (p2 - b^)u) 01. (3.3) 
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The forward scattering problem, explored in the next Section, is the problem of solving equations p.2p and 

on r. 



3.4. Boundary Conditions. 

The boundary conditions that we assume on the soluton u are those exhibited by the known soliton 
solutions [2F of Q35. If we define 



= 1 V(Pr - a){pr 



^) iTM ' ■= i + p{k) ^'-'^^ 

then a one-soliton solution to Q3 (which depends on all N lattice variables) is given by 

u ^AT{a, b) S{a, 6) + S J"(a, -h) S{a, -h) + 6 J"(-a, h) S{~a, b) + V T{-a, -&) 5(-a, -6), (3.5) 
where the four constants are restrained by 

A'Dia + bf -'Bt(a-bf = — . (3.6) 

loao 

The corresponding dual function XL is 

IX = (a + 6)yiJ'(a, fe) V{a)V{b) + (a - 6)S J"(a, -&) V{a)V{-b) 

~ {a-b)QF{-a,b)V{~a)V{b) - (a + &)D J-(-a, -&) V^(-a)y(-6). 

Using this as a guide, and bearing in mind the choice a > & > 0, the boundary conditions that we assume on 
the solution u are 

u QF{-a,b) + T) F{~a,-b) as i -00 (3.7a) 
u ^Xo-AFia^b) + 'Ki'i>F{a,-b) as i -J- +00, (3.7b) 
where %o and %i are constants and the plane-wave factors J- satisfy 

V(p(i) -a)(p(i) -5)7 
The corresponding boundary conditions for U are 

U {a-b)eT{~a,b)-{a + b)T)F{-a,-b) as i -00 (3.9a) 

Ur^Xo{a + b)AF{a,b) + Xi{a-b)'BF{a,-b) as i ^ +00. (3.9b) 

4. Forward Scattering of 0i 

We are now in a position to carry out the forward scattering problem. We firstly consider the forward 
scattering problem for 0i, which is determining the solution of p.2p . given the above boundary conditions. 
With this choice of boundary conditions, by explicit calculation it follows that the function Q, defined by 

u 

has the following asymptotic behaviour: 

fi^p + p as i ±00. (4.2) 

This asymptotic result shows that the object behaves like a difference of Hl-type soliton solutions as 
i — > ±00 (see e.g. [ID]). This is perhaps not surprising as in [25] and [24] Miura-type relations between 
soliton solutions of HI and QSs were found, which take precisely the form of the quantity fi. This object is 
exactly the term appearing in p.2p . and we may therefore rewrite this equation as 

(p' - C')^li - (p + P + t^) 01 + (p' - C')^ = 0. (4.3) 
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Definition 4.1. Given an initial condition u — u{i) along the staircase T , the potential v = v{i) is defined 
to be 

+ =u:=l]-p-p= f j-p-p, (4.4) 

where U is determined by (|2.4p . 

With this definition the potential v vanishes at either end of the staircase F. 
Definition 4.2. The Jost solutions ip^tp to (|4.3p are defined by the boundary conditions 



i-i 



r=0 
i-1 



'^(^^ - n ( ^I^pI j " as z ^ -oo (4.5a) 



^(^;C)-n(^|^j' asz^-oo, (4.5b) 
and the Jost solutions "01 V' equation (14. 3p are defined by the boundary conditions 



'1 



r=0 
i-1 



^(^; ~ n ( pj^j+^ j " as * ^ +00 (4.5c) 



^(«; C) ~ n (^ pj^j^^ j ' as z ^ (4.5d) 

Since equation (|4.3|) is invariant under the map C — s> — C, it follows by the definition of the boundary 
conditions for the Jost solutions and uniqueness of the boundary value problem [22] , that 

^{i-0 = v{i\-0, ^(«;C) = V'(«;-C). (4.6) 

Since the general solution to (|4.3p involves two linearly independent solutions we may write 

V' = A(^ + B<^, ?/> = A<^ + B(^ (4.7) 
where A and B are independent of i and k{C,) = A(— C) and B{() = B(— C). 
Proposition 4.1. If is purely imaginary then 

|A(C)P = l + |B(C)p. (4.8) 

Proof Firstly given any two solutions x{i) and y{i) of (|4.3p . by eliminating the potential term one can show 
that the Wronskian 

Wix,y):={p^-e)Hxy-xy) (4.9) 



is independent of i. Furthermore if C is purely imaginary then equation (|4.3p is purely real and thus ip* and 
^* (the complex conjugates of (p and ^) are also solutions of this equation. By comparing the boundary 
conditions for ip* , tp and ip* , ip, by the uniqueness of the boundary value problem we have 

^*(z;C) = ^(z;C), rir,c) = ki;C)- 

By taking the complex conjugate of (14. 7p we then have A(C) = A*(C*) and B(C) = B*(C*). Now due to the 
linearity and anti-symmetry of the Wronskian we have 

W{iP,iIj*) = W{kip* +B(p,k*ip + B*(f*) = (|A|2 ~ \B\^)W{ip*,ip), 

and since the Wronskian is independent of i, these may be evaluated at the relevant boundaries which gives 
Wiip, Ip*) = Wi^p^ip) = 2C. This proves gH). □ 
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4.1. Analyticity and Asymptoticity Properties of the Jest Solutions. We now determine asymp- 
toticity properties of the Jost solutions as functions of the discrete independent variable i, and analyticity 
and asymptoticity properties of the Jost solutions as functions of the spectral parameter 



Definition 4.3. The functions A, A and T,T are defined by 



^(z;C)=A(*;C)n 



^(z;C) = T(z;C)n 



pM - C 



r=0 
i-1 



^(*;C)=A(*;C)n 



, ^(*;C) = t(*;C)n 



L^n VpW + C 

p(r) + C 



r=0 
i-1 



,L^o VpW-C 



Proposition 4.2. For C 7^ the functions A and T satisfy the following summation equations: 



1-n 



T(*;C) = 1 + 



-. +00 

- F 



in 



P{r) - C 

pW + C 

pW - C 
■i VpM + C 



vil)A{l;C) 
v{l)T{l;C). 



(4.10a) 
(4.10b) 

(4.11) 
(4.12) 



Proof. Equation (|4.3p for A becomes 

(p + C)^ - (p + p) A + (p - C) A = UA 
which in terms of i may be written as 

p(z + 1) [a(z + 2; C) - A(i + 1; 0] - p(i) [A(i + !;()- A(z; C)] + c[a(« + 2; C) - A(i; C)] = u(i + l)A(i + 1; C), 
and after summing from I =^ —oo to I = i — 1, and using A—)- 1 asz-T'— oo, we have 

i 

A(z + 1; C) [p(*) + C] - A(z; C) [p(*) - C] = 2C + ^ ^(')^(^; 0- 

which gives 

A(* + 1; Os{i + 1) - A(i; C)s(z) = [^(i + 1) - s(z)] +- [5(2 + 1) - s(i)] ^ ^(OA(Z; C) 



We now multiply this equation by the summing factor s{i) 11^=0 ( p(r)-^ ) ; 



and then we sum from j = io < i ^ 1 to j = ?' — 1, obtaining 



A(z; C)s(z) - A(z„; C)s(zo) = [s(z) - s(z„)] +- ^ [^(j + 1) " ^(j)] E ^(0A(^; 0- 

We now let io — >■ — cxd and assume that s{i) as i ^ —00. By changing the order of summation the double 
sum can be rewritten as 



E 1 

j— — oo I— — 00 



E E [<? + l)-5W]^(0A(?;C)= E ^(0A(?;C)E[^O' + l)-^W 

'=-OD j = l 

E [s{i) ~ s{l)]v{l)A{l;C) 



and so the summation equation becomes 



A(^;C)s(^)=s(^) + - E [s{^) - sil)]vm{l; () 



l — ~oo 



A(*;C) = l + ^ E 



I — — 00 



1-n 



p(^) - c 

p(r) + C 



t;(OA(/;C), 
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which is equation (14. lip . Equation (j4.12p follows in a similar manner by summing to z = +00 and using 
T — s- 1 as i +00. □ 

Proposition 4.3. At = the functions A and T satisfy the following summation equations: 



A(z;0) = l+ J2 



/ — — OC 



T(z;0) = l+ J2 



=i+l 



l-l ^ 



(4.13) 
(4.14) 



Proof. At C = 0, equation (|4.3p for A becomes 

p{i + 1) A{i + 2; 0) - A{i + 1; 0) - p{i) A{i + 1; 0) - A{i; 0) = v{i + l)A{i + 1; 0), 
which after summing from I ^ —00 to I ^ i — 1 gives 

i 

p(*)[a(* + 1;0)-A(z;0)] = ^(0A(/;0). 

/— — 00 

After dividing through by p(i) and summing again from j = — 00 to j = i — I we obtain 



A(z;0) = l+ J2 E 
J— —00 /— — 00 

and by changing the order of summation we have 



1 



A(z;0) = l+ J2 



l=-c 



t-1 ^ 



v{l)Ail;0), 



v{l)A{l;0) 



□ 



which is equation (|4.1ip . The result (|4.12l) follows in a similar manner. 
Proposition 4.4. For C 7^ the summation equations (j4.1ip and (|4.12p have the Neumann series solutions 



+00 



A(^;C) = E^^ T(,c) = E 

k=0 ^ fc=0 ^ 



Jfc(»;C) 



where 



1 

Ho = l, i?fc+i(*;C)- 2 E 

/ — — 00 

+00 

Jo = l, Jfe+i(z;C) = ^E 



in 

l-l 

1-n 



p(0 - C 

p(r) + C 

p(0 - c 



^^(0^rfe(?;C), 



^;(Z)Jfe(/;C). 



(4.15) 

(4.16) 
(4.17) 



Proof. Inserting this series expression for A into the summation equation (|4.1ip gives 



A(.;C) = i + ^ E 



in 



r=l 

+00 ^ / ^ i-1 

= 1 + E7ITT 2 E 
^ -, gfc+i(»;C) 



p(0 - C 
pW + C 

in 



' +00 



^(0 E 



vil)Hk{l;C) 



k=0 



^Hk{i;0 



k=0 



as required. The proof for T follows in a similar fashion. 



□ 
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Proposition 4.5. At = the summation equations (|4.13p and (|4.14p have the Neumann series solutions 



+ 00 



+ 00 



fc=0 



k=0 



where 



/ — — CX3 
+ CXD 



Jo° = l, J,°+l«= E 
i=i+l 



P(j) 



^ p(j) 



3=1 



Proof. Inserting this series expression for A into the summation equation (|4.13p gives 

/ + CXD \ 



i-l 

A(z;0) = l+ ^ 

^ — — OC 


§^ 


+ 00 / 


i-l 


i-1 




E 

— — oo 


E 

j=i 


+c>o 






= i+E^^ 


















= E^^°« 







p(j) 



(o|E^fe(o 

>fc=0 / 



k=0 



as required. The proof for T is similar. 
Theorem 4.6. Assume that 



+ 00 



^ K*)|(l + |.|)<oo, 



(4.18) 

(4.19) 
(4.20) 



□ 



(4.21) 
(4.22) 



and that p(r) > for all r e I. Let TZ^ denote the half -plane 

Tl+:^{C : Re(C)>0}. 
Then for C e 

|A(i; C) - 1| < C'l for C ^ (4.23a) 

|A(i;C) - 1| < C2(l + max{0,i}) (4.23b) 

|T(i;C)-l| <C3 for Ct^O (4.23c) 

|T(i; C) - 1| < Ci{l + max{0, -i}) (4.23d) 

where Ci — > C4 are constants. For all C, G 7?.+ i/ie series solutions for A and T converge absolutely in i, and 
uniformly ifC,^0. For each i, A and T are continuous functions of in TZ^ , and analytic functions of in 
the interior of this half-plane. 

Proof. The proof of this Theorem is obtained by showing absohite and uniform convergence of the Neumann 
series representation of the Jost solutions. The complete details, reminiscent of the analysis given in |13) for 
the continuous case, are given in the Appendix. The estimates obtained agree with those obtained in [5], [32] 
and [in] for the a similar spectral problem. □ 



Corollary 4.7. For C £ 7^+ the functions A and T have the following asymptotic behaviour: 

'1~ 



A(z;C) = l + 
T(z;C) = l + 



as Id — 00 
as Id ^ 00. 



(4.24a) 
(4.24b) 
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Proof. From the series solution of A we have 



+00 



A(*;C) = 1 + E 



gfc(»;C) 



fc=i 

For C G TZ^, C 7^ however, from the Appendix we have that \Hk\ < K for some constant K, and thus 
Hk = 0{1) as \C,\ — > 00, for all k > 1. This proves the result, and a similar argument works for the series 
solution for T. □ 

Corollary 4.8. Theorem \4.6\ and Corollary \4.'T\ hold for the functions A and T for C, in the half-plane 

n- :={C:Re(C)<0}. 

Proof This follows from the fact that k{i; () = A{i; -Q and f{i; C) = T{i; -(). □ 

4.2. Analyticity and Asymptoticity Properties of A and B. We now look at analyticity and asymp- 
toticity properties of A = A(C) and B — B{(), which are defined by equation (|4.7p . These are related to the 
reflection coefHcient R and transmission coefficient T by 



A 



Proposition 4.9. A and B have the following properties: 

- k is analytic in the interior of TZ^ and continuous in TZ^ , except possibly at — 

- B is continuous on the imaginary Q-axis, except possibly at ( = 0. 

Proof. Taking the Wronskian of ip ^ Aip + Bip we have 

A(C) = ^Wi^P, ip) = ^ ((p(*) + C)A(* + C)T(*; C) - (p(^) - C)A(*; C)T(^ + 1; C)) 



2C 



p(»)-C 

2C 



(a(*;C)T(z + 1;C)-A(j + 1;C)T(*;C)) 



Since A and T are continuous in TZ^ and analytic in the interior of this region, A also has this property, 
except possibly at ^ = 0. The expression for B however is only valid on the intersection of 7^+ and TZ~ , 
i.e. the imaginary ^-axis. Since A and T are continuous here, B also has this property, except possibly at 

C = o. " □ 



Proposition 4.10. For C 7^ the functions A and B can be expressed as 

+00 



, +00 

^(0 = ^ E 



pM - C 



7-1 

n Vp(r) + C 



vil)Til;C). 



(4.25) 
(4.26) 



Proof. The summation equation (I4.12p for T may be written as 



T(*;C)= 1 + ^ E ^(OT(?;C) +n 



p(r) - c 

V{r) + C 



^ +00 
2C ^ 



n 

.r=0 



p(0 - c 
p(^) + c 



Taking the limit i — >■ — cxd and comparing this with 



v{l)T(l-C) . (4.27) 



T(z;C)'-A(C)+B(C)nf^|4T7 



as I ^ —00 



gives the desired result. 



□ 
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Proposition 4.11. For C G 7^"*" we have 



A(C) = l + 0(-) as ICI^oo 



and for C on the imaginary axis we have 

B(C) = o 



as Id oo. 



(4.28) 



(4.29) 



Proof. Inserting the asymptotic behaviour (|4.24p of T into the expression (|4.25l) for A gives 

+ 00 



/ — — oo 



and since 



+ 00 



E "(0 



+ 00 

< E H0\<^ 

1— — 00 



this proves (|4.28p . Performing the same task for the expression (|4.26p for B gives 



7-1 



B(o = ^ E 



n 



P(0 + C 



l + O 



and since for purely imaginary ^ we have 



+ 00 

E 



■;-i 

n 



p(0 - c 

p(r) + C 



vil) 



+ 00 

< E 1^(01 < 



this proves 



□ 



Theorem 4.12. T/ie function A /las a finite number of bounded isolated zeroes {C/j, /c = l,...,Af} m t/ie 
interior of 72."*" , and moreover every Qk is purely real and satisfies Cfe < Pr /o?' parameters pr existing along 
r. At each zero of k we have '4>{i',Ck) — ^{Ck)^{i'Xk), md 



+ 00 

E 



'fii-XkMt + i-Xk)] _ A'(a) 



where denotes the derivative of A with respect to C.- 



HCk) 



(4.30) 



Proof. Since A ^ 1 as |C| — oo it foUows that there exists some constant Co such that < Co for every k. 
Since 2(k{() = WiiJjjLp) it foUows that for every fc, (p{i;(k) and V'(i;Cfc) s-re hnearly dependent, so we may 
write tl){i]C,k) = bk(p{i',Ck) for some constant bk. This imphes that 



i-l 



T(*;a)^&fcn 
and so by equation (|4.27p we have 



p(0 - g 



as z — !■ — oo, 



/ — — oo 



^ +00 



i-l 

n 



p(0 - Ck 

P{r) + Ck 



vil)T{l;Ck). 



Thus bk — B(Cfc) for every k. Now consider the scattering problem (|4.3p for at ^ = Cfe: 

(p' - Cfe)^ ^(z + 2; Cfe) - (p + p + u) ifii + 1; a) + (p' - Cfe)^ ^(i; a) = 0. 



A DISCRETE INVERSE SCATTERING TRANSFORM FOR Q3s 

For every (^k in the interior of 7?.+ we have 



13 



4, VpW - Cfc 



r=0 



as « 



^(z;a)~B(Cfc)n 



Pir)-Ck Y 

p(0 + a 



— ?> as i — ?► +CXD, 



and thus (p{i; (k) is summable over all i. If we multiply the scattering problem for ip{i; Cfc) by (p*{i + 1; Q), 
sum over all i and define 



s(i;Cfe) :=<^(i + i;Cfc)¥'*(«;Cfc) + 'P*(« + i;Cfc)<^(«;Cfc) e 



then we have 



+ 00 



+ 00 



E (p«'-^')^^(*;Cfc)= E [pW+p(i + i) + «('; + i)]i^(* + i;Cfc)P, 

i— — oo i— — oo 

which implies that Cfc must be real and that < Cfc ^ p(*) for every i. Thus Cfc must be less than every 
parameter pr through which p(i) cycles, which proves the given statement. 

Now A has isolated zeroes along the positive real C axis which are all bounded. The only way that there 
could be an infinite number of these zeroes is if they formed a limiting sequence which accumulated at C = 0. 
We will show that this is not possible. Suppose that such a sequence {Cfc} of zeroes exists: limfc_j.oo Cfc = 0- 
Then at each Cfc we have 



B(Cfc) 



and so 



lim |B(Cfc) - B(0)| = lim 



V'(i;Cfc) 



= 



Cfc) 

since the Jost solutions are continuous at C = 0. At C = however we have 1^9(1; 0) = ip{i; 0) and so 



A(0) +B(0) 



V'(»;0) 

^(z;0) 



lim 

k^oc 



V'(»;Cfc) 
<^(«;Cfc) 



lim B(Cfc) - B(0) 



fc— )-oo 



which implies that A(0) — 0, which in turn contradicts Proposition 23] Thus A has only a finite number of 
zeroes in TZ'^ . 

Finally to prove (|4.30p we define the following two useful functions: 
We differentiate 2CA = W{iIj, (p) to obtain 



2CfcA'(Cfc) = W{ij{i; Ck),v'{^■, Cfc)) + W{i>{i; Ck),^'{i; Cfc)) = HCk)W^{i; Cfc) 



B(Cfc 



-W^{i;Ck). (4.31) 



Now consider the difference of two equations: firstly the derivative of the scattering problem (|4.3|) for ip 
multiplied by 1^9(1+!; C), and secondly the (un-differentiated) scattering problem for ip multiplied by ip'{i+l; C). 
This gives 

M/r^i n w (■ t\ .r¥'(*;C)^(* + i;C) , v^C* + i;C)¥'(« + 2;C)" 

which may be summed to give 

i-l 

W'^(*;C) = C E 



(p(z)2-C2)5 (p(, + l)2_^2) 

^G;C)^a + i;C) , p{l + i-X)v{l + 'i-X) 



(pa + i)-c^ 



One can then perform the same task with ^, only instead this time summing from i to +00, to obtain 



+ 00 r 



w^v(*;C) = -cE 



i=i 



i;{lXm + ^-X) . V^(/ + i;C)V'(^ + 2;C) 



ip{ir-ey. 
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Now set C = Cfe and rewrite (k) ~ ^{Ck)'p{i', Cfe)- The expression (|4.3ip then becomes 



+00 



I — — CO 

+00 

= 2aB(a) J2 



ipil;CkMl + l;Ck) , </j(/ + l;Cfe)v(' + 2;a) 



I — — 00 _ 



(p(O^-Cl)' 

y)(/;Cfc)ya + l;a) 



(pa + i)-ci)^ 



which gives (|4.30p . 



□ 



The sum in equation (j4.30l) is of fundamental importance to the scattering problem, and as such we make 
the foUowing definition. 



Definition 4.4. The square eigenfunction $ is defined to be 

(pW'-C^)^ 

and we define the normalisation constants {c/,., k = l,...,Af} to be 

+00 



Cfc 



B(a.) 



-1 



A'(a) 

From this definition it is clear that we have the fohowing resuh 
Ttieorem 4.13. If the normalisation constants Cfe are all finite, then every zero A in TZ^ is simple. 
One possible way of ensuring that all normalisation constants are finite is by using the result 



(4.32) 



(4.33) 



+ 00 



ip{i-Xk)^ii + i;Cfc) 



+ 00 



= E [?i^)+v{^ + l) + vi^ + ^)M^■Xk)\'■ 



Thus if we impose, for example, that [p(z) + p{i + 1) + v{i + 1)] > for all i, then every will be finite. We 
do not dwell on this further, but assume henceforth that all normalisation constants are finite. 

This concludes the forward scattering of (j>i. The main results of this section are the construction of the 
Jost solutions A{i; (),T{i; the spectral functions A(C),B(C), and the knowledge of their analyticity and 
asymptoticity properties. 



5. "Time" Evolution of the Scattering Data 

We now consider how the spectral functions A, B and the normalisation constants depend on the N 
lattice variables. As an analogy to the continuous theory, this is the calculation of the "time" dependence of 
these functions, with respect to the arbitrary number of discrete "time" variables. 

From the Lax pair the equation governing the evolution of the Jost solutions (obtained by eliminating the 
second component) in any one particular lattice direction with variable Uk and parameter pk is given by 



(5.1) 



where u denotes an iteration of u in the n^-direction. From the boundary conditions (|3.7p for u, for all 
fc = 1, N we have 



U 



2pk as i — > ±cx). 



(5.2) 



In other words at both ends of the staircase the dependence of the Jost solutions on the lattice variable Uk 
is governed by 

{pI - eY^^- ^Pk?+ (pI -e)h^o- (5.3) 
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5.1. "Staircase Directions" (e I). 



Recall that I denotes the set of lattice directions in which the staircase of initial values V iterates. Due to 
the way it was defined in Section [31 the set of initial values along T defines a well-posed initial value problem 
[B] for all points in each of the "staircase directions" indexed by I. For example given the staircase of initial 
values in Figure [Sal all points in this plane are uniquely defined by (|l.ip . Likewise in Figure [5bl given this 
staircase of initial values, all points in this three-dimensional space are uniquely defined by (jl.ip . 




(a) All points are uniquely . s are unique y ^-^^ Points in the plane 

determined by T determined by i uniquely determined but 

those in the vertical direction 
are not 



Figure 5 



For all "staircase directions", in order that the Jost solutions be consistent with each of the equations 
5.31) . we redefine their boundary conditions (|4.5al) (|4.5cP as follows: 




LP 



n ( ^ j as z ^ (5.4a) 



Pr ~ C 
Pr + Q 

Pr + C 
Pr-C 



as i +O0. (5.4b) 

\ — [ / 

When one restricts these boundary conditions to F these agree with the previous boundary conditions (|4.5ap 
(I4.5cp . Now for each "staircase direction", (p and (p are linearly independent functions of equation (|5.ip . so 
the functions A and B, and the normalisation constants c, are independent of all "staircase variables" Uk, 
kel. 

5.2. "Non-Staircase Directions" T). 

We now consider how the spectral functions A and B and the normalization constants c evolve as we move 
in any of the "non-staircase" lattice directions (^ I). We denote the collection of these "non-staircase" 
direction by JT", so that IVJ J \s the entire N-dimensional lattice. For example in Figure [Sc] the two lattice 
directions in the lower plane are elements of I, but the orthogonal vertical direction is an element of J . In 
order that the evolution of (jl.ip be well-defined in each of these directions we make the following assumption 
on the boundary conditions of u at one end of the staircase: 

In the limit i — ^ — oo, the boundary conditions (|3.7ap for u, and (j3.9ap for 11, are assumed to remain 
unchanged under any shift in a "non-staircase" direction J). 

This assuption gives sufficient data for the solution u to be determined uniquely at all lattice points using 
(|l.ip . This is most evident when we consider a line of initial values Fq as shown in Figure [H In this case I 
is the horizontal direction, and J is the set of all such orthogonal vertical directions. 

The assumption of the invariance of the boundary conditions of u and U (as i — oo) gives sufficient data 
for the solution to be iterated to all points in the lattice. The boundary conditions (|5.4ap (|5.4bp however are 
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T3 



a 
o 



o -a 



o 

u 

O 



o 



-oo 



Figure 6. Boundary conditions giving sufficient data to iterate the solution in the "non- 
staircase" vertical direction 

assumed to hold independently of these "non-staircase variables" Ug . Since the evolution of the Jost solutions 
in any one "non-staircase" direction is given by (|5.ip . our boundary conditions are inconsistent with (j5.3|) . 
To circumvent this we make the following definition. 

Definition 5.1. Let J denote the collection of "non-staircase directions". The N- dimensional Jost solutions 
(pW^ipW and il)'^^\'ip'^^\ which are solutions to all N equations (j5.1l) . are defined to he 



n 



Ps+C 



sej 

(N) TT (Ps~C 



Ps~C 
Ps+C 



seJ 



Ps+C 



(N) J. TT /Ps+C 



Ps-C 



(5.5a) 



(5.5b) 



We then have the following result about the dependence of the spectral functions A and B and the normal- 
isation constants c on each of the "non-staircase variables" Ug, s ^ J . 

Theorem 5.1. 



The function A is independent of all lattice variables.- 

A(n,;C) -A(C) 

The dependence of the functions B and c on n^, s ^ J , are given by 

'Ps+C " 



B(n,;C)-B(C)n 



c(n,;C) = c(C)n 



Ps-C 

Ps+C 
Ps-C 



(5.6) 

(5.7) 
(5.8) 



Proof. Consider the "non-staircase" ris-direction, where s & J . The evolution equation for the N-dimensional 
Jost solutions in this direction is 

{pI - C^)^^^^' - (^^^) + (pI - C^)^^^"'^ = 0, 

and by their boundary conditions the functions ip^^^ and ip^-^"^ are linearly independent solutions of this 
equation. We may therefore write 
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where Ci and C2 are independent of n^, but may depend on all other lattice variables. This is equivalent to 



sej 

and by comparing this with 

>r + C 



rGl 

it follows that A is independent of n^, while the dependence of B on Ug is through the plane- wave factor 
■ Repeating this argument for all s G gives the desired results for A and B. The result for c 
follows from its definition (|4.33p . □ 



6. Inverse Problem for (j)i 

We now consider the inverse problem for 0i, that is the construction of the Jost solutions as functions of 
all N lattice variables. In doing so we alter the notation of all eigenfunctions by writing 

and similarly for the other functions. Here (p{C) is understood to depend on all lattice variables ni,...,njv, 
however for the inverse problem it is convenient to suppress this dependence in the notation. 
Consider the equation (j4.7l) . which we rewrite as 

T(C) 



A(C) 

where the reflection coefficient R is given by 

B(C) 



-A(C) = i?(C)A(C)p(C), (6.1) 
RiC) = 



A(C) 

and the plane- wave factors p are defined by 

This defines a jump condition between two sectionally meromorphic functions along the contour Re(C) = 0. 
To 

The functions — and A are analytic in the interior of the regions TZ^ and TZ respectively, and both are 

continuous along Re(C) = 0. Given the jump condition and their boundary conditions, the question of 
determining a function which is equal to these in their respective half-planes is a Riemann-Hilbert problem. 
The method of solving such a problem is well-known (see e.g. Section 4): Consider the singular integral 



7^/ \ Pi<^) dcr^ — , ' .^ da- — \ da (6.2) 

T 

where Q G TZ^ . Here the contour of integration is the imaginary cr-axis. Since the function — has M simple 
poles in 7?.+ and has the boundary behaviour -j- ^ 1 as |<^| — 00, one may use the residue theorem to calculate 

1 /•+*°° T(a) ^ 1 _ _ T(a 

' ^ 2 ^ A'f 



2^1 J k{a){cT + Q) 2 ^A'(a)(C + a) 



By then using the fact that 
this can written as 



T(a) = B(a)A(a)p(a) 



27rj 



M 



1 /■+'°° T(a) ^ 1 ,^c,A(C 



A(a)(. + C)'"^2-l^(CTa)^(^'=) 
can determine that 

/ ^da = -- + A(-C) = --+A(C), 
J-ioo cr + C 2 2 



Secondly since A is analytic in TZ one can determine that 

27ri 
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and thus the singular integral equation becomes 

Given the Cfc, the normalisations constants c^. and the reflection coeflicient R, one can use this equation to 
determine the Jost solution A as a function of C. and all N lattice variables. Remarkably all of the dependence 
on the lattice variables is contained in the plane-wave factors p. If instead we started with the relation 

A(C) = A(C)t(C) - B(-C)T(C)p(-C), (6.4) 

which is consistent with (|4.7p . then by following the same procedure as above one finds that T(C) is determined 
by solving the singular integral equation 

where 

7. Inverse Problem for (/)2 

Now that we have constructed the first component of the eigenfunction <p of (I2.3p . we use the first com- 
ponent of the Lax equations to determine 02 : 

ie-b')MC) = (pI-O^umO - (Pku- ipl-b')u)MO- (7.1) 

Let if'^'^ and '02^'' be the corresponding second components for the Jost solutions (^(^) and ?/;(^) respectively, 
and define the functions A2 and T2 by 



N . s^n^ N , 



Then A2 is given by 



and T2 is given by 



Pr - CJ \Pr + C 



(C' - b^) A2(C) = {Pk + C) W A(C) - (Pku - ipl ~b^)u) A(C), (7.2) 



- ^'')T2(C) = (Pfc - C)UT(C) - (Pku- (pI - 52)z.)T(C). (7.3) 
We see that A2(C) and T2(C) are analytic in 7?.+ , except for a simple pole at C = Thus A2(C) = A2(— C) and 
T2(C) = T"2(— C), which are the second components for the Jost solutions A and T respectively, are analytic 
in TZ~ except for a simple pole at C = ~b. Furthermore all of these functions are continuous on the imaginary 
C-axis. Now since the two eigenfunctions 



are linearly independent solutions of the Lax equations (12.31) . we may write 

V''^^ )=A(C)<^W+B(C)<p(^), 
whose second component may be written in terms of A2 and T2 as 

^-A2(C)=i?(C)A2(C)p(C). (7.4) 

Note that these are the same functions A and B as those that appear in the integral equation (j6.3p for A(C). 
As in the inverse problem for 0i, equation (j7.4l) becomes the jump condition for a Riemann-Hilbert problem, 
and as such we look at the singular integral 

1 r*'- ^ 1 f T.(.) ^_>_r-AM ^ 
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The difference in this case is that T2 and A2 have simple poles at +b and —b respectively, and both of these 
function are 0(C~^) as |C| ^ 00. Let us first consider the integral involving T2. By using the residue 
theorem one has 



1 



k{a){a + 



mic+b) 



k=l 



(c + a) 



From the Lax equations however, by eliminating the first component (j)i one can show that at C ^ the 
second-order linear equation for (/)2 in the nfe-dircction drastically simplifies to 



Pk{u~u) 
U 



02 + (j3fc-a')3(/)2=O, 



(7.7) 



which we identify as equation (jS.ip aX Q = a. Since this holds for every lattice direction we may write 



for some constants a and j3. In terms of T2 and A2 this implies that 



N 



=a k{a) T{a,b) + l3T{a) T{-a,h), 
where the plane- wave factors J- are defined by 

N 



N 



-/?T(a)n 



— n 

{Pr - a) {pr + &) V 



(7.1 



J-(a,6) 



{pr + a){Pr + b) 
rJl \{Pr - a){Pr - b) 

In order to determine a and (3 we consider equation (j7.8p in the limit i — > ±00. Firstly we have the asymptotic 
behaviour 

A(a) ^ A(a) + Cop{—a) as i ^ +00 
T(a) ^ A(a) + cip(a) as z — > —00, 

for some constants Cq and ci. Since a > both of these plane- wave factors are exponentially small. Thus 
using the boundary conditions of u and U in (j7.3p . as i — > — cxd we have 

(C - 6)T2(C)|^^, - ^ [{Pk -b)U- {P,u- {pI -b^)u ) 
- {a-b)k{b)QT{-a,b), 
which implies that j3k{a) = {a — b) Ck{b). By then taking i +00 we have 

1 

IC^fc ' 2b 



(C-fc)T2(C)| 



{pk-b)U-iPkU-{pl-b^)u) 

{a + b)XoAJ'{a,b), 

which gives ak{a) = — (a + b)3Co-A. Therefore the integral in equation (j7.6p becomes 



1 

27ri 



T2(a) 



k{a)ia + 



da = {a + b 

M 

-E 



3C„ 



k{a)k{b) 



AF{a, b) 



(Ma). 

\C + b 



{a-b)QT{-a,b) 



T(a) 



A(a)(C + 6) 



CfcA2(Cfc' 

(c + a) 



■p(Cfc)- 



The integral involving A in (|7.5I) can be evaluated to be 



1 



- lim 

Q — b a^-b 



A2{a){a + b) 



(7.9) 



and by similar reasoning and using the fact that A(— a) = A(a) etc., we find 



lim 



A2{a){a + b) = -{a + b) 



Xik{b) \ 
A(a) J 



•BTia, -6)A(a) + (a + 6)1) J"(-a, -b) 



T(a) 
A(a) 
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Thus the integral (17.91) becomes 



1 



Using these in (j7.5p then gives the fohowing closed- form singular integral equation for A2: 
A natural ansatz for this is 



A2(C) = (a + 6) ) ^Ha, &)A(a)6(C) + («-&) (^-^ j 55-F(a, -6)A(a)e-b(C) 

- (a - fe)e^(-a, b) (^I^) 6(C) - (« + -&) ^-^(C) (7.13) 

where by equation (|7.12p the functions £,±b{0 are calculated by solving 

Note that the ingredients in these equations are the scattering data from the forward scattering of (pi, and 
(|7.14p differs from (|6.3p only in the source term. 

Remark 7.1. If we define the quantity 

S{a,b):=^-Ua) (7.15) 
a + 

then by comparing (|7.14p with the integral equation considered in the direct linearization approach [26] we see 
that this object is in fact a solution (containing solitons and radiation) of the NQC equation. 

8. Reconstruction of the Solution of QS^ 

We now show how one can recontruct the solution m as a function of all N lattice variables. Consider 
equation (|7.2p . which holds for any n^, fc = 1, A^. Firstly by dividing through by {pk + C) and taking the 
limit |(^| — > 00 we have 

11= lim [CA2(C)]. (8.1) 
Using equations (I7.13P and (|7.14p this may be expressed as 

" = + {wim) + - ^) -^)A(«)^(-^) 

- (a - 6)e^(-a, b) V{b) -{a + b)'DT{-a, -b) F(-&), (8.2) 



(8.3) 



where V{±b) is given by 

-1 /•-t-icxD 



V{±b) ^l-Y. Cfc^±fa(a)p(a-) - TT- / R{a)^±b{a)p{a) da. (8.4) 
fc=i 



In fact it turns out that A and V are the same object. Equation ()7.2p could then in principle be summed to 
find u. There is a way however to obtain a closed-form expression for u rather than its derivative. To do this 
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we first note that in the Lax equations (|2.3p . one is free to interchange the roles of the parameters a and h. 
In other words the N equations 




PkU 



(8.5) 



are also N Lax equations for QSd". Note however that we are not sqapping a and h in the functions u and 11, 
nor in the initial conditions, nor in the boundary conditions. We are simply repeating the entire 1ST with 
the new Lax equations ()8.5|) in place of ()2.3|) . which is permissible due to the symmetric dependence of QS^ 
on a and h. There are however some important remarks to be made. Firstly the forward scattering problem 
(j4.3l) for is independent of a and 6, and thus interchanging a and b in the Lax equations will not change the 
Jost solutions A and T, nor the scattering data. The second component (j32 however will now be calculated 
by swapping a and b in the Lax equations (|7.2p and (I7.3p . and re-deriving the integral equation (|7.14p . To 
make a clear distinction between the original quantities A2,T2 and their new counterparts, we rewrite the 
original functions as 



and then denote the new functions, which are obtained by swapping a and b in (j7.2p and (|7.3|) . by Aj"' and 
Tj"'. By repeating the analysis of the previous section we find 

^^'^^^ ^ + ^) (i(^) + - ^) (^) -'^ im) ^'^^^ 

-{a- b)QT{-a, fe)A(6)f_<,(C) - (a + b)'DF{-a, -b) e-a(C)- (8.6) 



7) 



A(fe) 

Then finally given the quantities K^'^ and Aj""* one can combine the two versions of (|7.2p , namely 

(C^ -62)Af (C) = (p, +C)UA(C) - [PkU- [pi b')u)A{C) 
- a') h^^^ (C) = {Pk + C) U A(C) - {PkU [pI -a')u ) A(C) 
to express the N-dimcnsional solution of Q3s as 

" (a2-62)A(C) 

Since however the solution is independent of ^, we may take it to be large, in which case the solution can be 
expressed as 

„,„ (cnAf'(o-f(o]y 

This solution is obtained by solving the integral equations (|6.3p . ()6.5p and (|7.14p for A, T and ^, using these 
objects in the expressions (|7.13p and (|8.6p for A^^^ and A*^°' respectively, and then taking the limit in (|8.8p . 

Remark 8.1. 

For the reflectionless case B = the Jost solutions exist at the points ^ = —a and ^ = —b, and we have 
Xo = A(a)A(6) and 3Ci — A(a)/A(fe). Comparing the expressions for u obtained by setting ^ = ±a, ±& in (|8.7p . 
and using the fact that T(C) = A(C)A(— we find 

u = {a+b)AT{a, b)ia{b) + {a-b)'BT{a, -b)ia{-b)-{a-b)QT{-a, b)C-a{b)-{a+b)VT{-a, -6)C-a(-6)- (8.9) 

Using the NQC variable S{a, b) defined by (|7.15p this becomes 

u ^AT{a, b) [1 - (a + b)S{a, b)] + S J"(a, -b) [l - (a - b)S{a, -b)] 
+ ej'(-a, b) [l + (a - b)S{-a, b)] + TiT{-a, ~b) [l + (a + 6)S'(-a, -6)] . (8.10) 
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This is precisely the form of the N-soliton solution of QSs obtained in |25j . The corresponding dual function 
U is given by 

U = (a + b)AJ'ia, b)A{a)A{b) + (a - 6)'B J"(a, -b)A{a)A{~b) 

- (a - b)ej-{-a, 6)A(-a)A(6) - (a + b)T>T{-a, -b)A{-a)A{-b), (8.11) 
which was also obtained in (25) . 

9. One-soliton Example 

Here we construct an explicit solution to QS^ for the special case where the reflection coefficient R is 

identically zero and the function A has exactly one zero in TZ^ at — k > Q. The constants 3Co and "Ki are 

set to ((J^^m^) and 
SLL LO J anu 

(|7.14p for We write the normalisation constant as ci 
into the plane-wave factor p. Equation (|7.14l) becomes 



^ ^ respectively. We are required to solve the singular integral equation 

2kc with c being a constant, which we absorb 



which has the solution 



1 



2k 



C±b \c+k 



C±b 



From 



the solution is then 
1 

u — AT{a, b) 



+ ej'i-a,b) 



{a-k){b-k) 
(a+k){b+k) 



p{k) 



l+p{k) 

(a+k)(b-k) 
\a-k)(b+k) 



{C-k){bTk) 
JX+k}(b±k)f}_) 

l+p{k) 



^T{a,~b) 



p{k) 



1 + p{k) 



(9.1) 



(9.2) 



(a-fc)(b+fc) 
(a+fe)(6-fc) 



p{k) 



l+p{k) 



{a+k)(b+k) 
(a-k)(b-k) 



pik) 



l+p{k) 



(9.3) 



where AV{a + bf - Se(a - bf = 



10. Conclusion 



In this chapter we have rigorously derived a discrete 1ST for the QS^ lattice equation. The initial-value 
space was given on a multidimensional staircase within an N-dimensional lattice, and we have given examples 
of how such staircases may be constructed. By incorporating the multidimensional consistency of the QS^ 
equation into the 1ST scheme we have shown how to solve the inverse problem for this equation for solutions 
depending on an N discrete independent variables. The assumptions made on the solution were that it be 
real and that the initial profile satisfy the summability condition ()4.21|) . which is a much weaker condition 
than that imposed in [9]. The solution to QS^ incorporating both solitons and radiation is found by solving 
the singular integral equations (|6.3p . ()6.5p and ()7.14p . and then taking the limit given in ()8.8p . This method 
of solution should of course also apply to the other ABS equations, and will form the basis of future research. 

11. Appendix 



11.1. Proof of Theorem 14.61 We first prove ()4.23a|) . For ^ G 7?.^, C 7^ the recursion relation (|4T6|) can 
be upper-bounded by 



i-l 



We then claim that 



where 



\Hk+i{i-x)\ < E HmHkii-x)\ 

1 — — 00 

Fit)'' 



\Hk{^-X)\ < 



k\ ' 



(11.1) 
(11.2) 
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Clearly this holds for fc = 0. To prove the inductive step we use summation by parts and the fact 

that F{i + 1) > F{i): 



/ — — oo 

= ^ E [F{l + l)-F{l)]F{ir 



1— — 00 

k+1 



< 



fc! 
fc! 
fc! 



fc! 



-.1-1 / fc-i 

- E [F{l + l)-Fil)]F{l + l)i Y,Fil + l)'^-'-^FilY 



1— — 00 

i-1 



r=0 



fc! 



E [F{l + l)-F{l)]F{ir, 



and so by examining the second and last lines we have 



1 — — 00 



(fc + 1)! 



which then shows that 



\Hk+i{^;0\ < 



F{^) 



k+l 



(fc + 1)!' 

Thus the estimate (|11.2p holds. The series solution for A can then be upper-bounded by 



|A(^>c)-ii<E^^<E^<(^^ 



exp 



F(+oo) 



ICI 



<Ci 



since (|4.2ip holds. Thus (I4.23al) is proved and so for any C 7^ the series solution for A converges absolutely 
and uniformly in i. We now prove ()4.23bp . To allow for C — give an alternative upper-bound for the 
summation equation (|4.1ip . One can easily verify that 



1-n 



r=l 



P{r) - C 

pW + C 



E 

3=1 



p(j) - c 

p(.?) + C 
pM - C 



J-1 

n 



i-1 / i-1 



< 



1 + E I E ip(j-) + ^1 



pW + C 

v{l)\\A{l;C) 



p(0 - C 
pW + C 

^;(0||A(^;C)| 



(11.3) 



/ — — OO — l 

For G TZ^ however we have 

m + CI > ip(j) " CI > I ip(j)i - ICI I > ip(j)i - ICI > ip(j)i 

and so we have 



|A(*;C)|<l + ^7 E (» -0 1^^(01 |A(^; 01, 



(11.4) 



where rj = max{ |p(r) ^\ : r € X}. Equation (|11.4[) is a majorant for both the summation equations (|4.1ip 
and ()4.13p . and thus may be used to estimate for all C S "T^"*"- Thus we have 



+ 00 



|A(z;C)|<E'?'^fe« 



k=0 



where 
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H*^ = 1, = ^ - \v{l)\Hl{l). 



We claim that 



where 



k\ 



(11.5) 



G{i,3)= E (^-OKOI- 



Clearly this holds for fc = 0. To prove the inductive step we use the recursion relation for H^, and again 
summation by parts and the properties of G: 



m+iii)\< {i-i)\v{i)\ 



1 — — 00 

2-1 



G{l,lf 
k\ 



[G{i,l + '^)-G{i,l)]G{i,lY 



1 — — 00 



^G(i,i)'=+i-i £ [G{i,l + l)'-G{i,l)']G{i,l + l) 



-oo 
i-1 



-G(i,i)'=+i-- J2 [G{i,l + l)-G{i,l)]G{i,l + l)[ J2G{i,l + l)>'-'-^G{i,iy 



L— — 00 

'^'-^ E [G{i,l + l)-G{i,l)]G{i,l)' 

I— — CO 



r=0 



and again by examining the second and last lines we have 



i-1 

- E [Git,l + l)-G{i,l)]G{i,l)'' < 



G{i,i 



|/C+1 



(fc + 1)! 



which implies 



G{i,i) 



k+l 



{k + iy. 



and completes the inductive step. Thus for all ( G 7^+ 



|A(»;C)-1| <E '^ u <vG{i,i)exp 7jG{i,i) 



k=l 



Let us first consider i <0. In this case we have 

|A(i;C) - 1| < ?7exp r]G{i,i] 



i-1 



< ryexp 



r?G(0,0)] 



r=— oo 
-1 
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for some constant Di, courtesy of (|4.21l) . To examine the case z > we consider the majorant (|11.4p 

i-l i-1 

m-X)\<i + v E H)KOI|A(/;C)l + »^ E KOI|Aa;C)l 

^ — — oo — oo 

-1 i-l 

<i + vJ2 H)Hi)\m-X)\ + ^v E 1^(01 |A(/;C)I 

^ — — oo — oo 

i-l 

<D2 + iiiYl Hi)\\Mi-X)\ 

1 — — 00 

for some constant D2, where we have used (|4.2ip and the fact that A(i; C) can be upper-bounded by a constant 
for J < 0. Write A(i; C,) = 1)2(1 + C), then the upper-bound for S becomes 

|2(*;c)i<i+^(y^) E {i + \i\)\v{i)\mc)\ 

i-l 

<i+v E {i + \i\)\v{i)\mo\- 

I— —00 

Making similar arguments to those presented above it follows that 

|S(z;C)|<expU E (1 + 1^1) KOI ) <cxph E (1 + 1^1) KOI ) 

\ I — — 00 / \ I — — 00 / 

which in tm^n implies that for i > 0, 

|A(*;C)I <C2(i + z) 

for some constant C2. Combining this with the result for i < proves (j4.23bp . Thus for each i the series 
solution for A converges absolutely and uniformly in ( for ^ G T?."*". Since the iterates Hk are continuous in 
this region and analytic in its interior, A also has this property. The results for T follow in a similar fashion. 
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